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The  potential  importance  of  energy-level  saturation  to  analytical  atomic 
fluorescence  spectrometry  has  prompted  a two-part  fundamental  investigation 
into  the  saturation  process.  The  first  part,  reported  here,  develops  a model 
for  saturation  which  is  applicable  to  both  resonance  and  non-resonance  fluores- 
cence of  atoms  in  flames  or  other  gaseous  media.  Provision  is  made  for  con- 
sidering both  collisional  excitation  and  quenching  of  atoms,  and  the  exciting 
source  can  be  chosen  to  be  of  any  desired  duration.  To  illustrate  the  utility 
of  the  model,  values  are  predicted  for  the  source  spectral  power  density 
necessary  to  saturate  several  analytically  important  atomic  transitions.  Ap- 
plication of  the  model  requires  only  a knowledge  of  the  kinetic  parameters  of 
the  energy  levels  being  probed,  including  Einstein  coefficients  and  quenching 
rates.  Assumptions  implicit  in  the  model  are  discussed  and  its  sphere  of 
valid  application  considered.  Experimental  verification  of  the  model  will  be 
reported  in  part  II. 


A number  of  investigators  have  explored  the  application  of  tunable  lasers 
as  excitation  sources  in  atomic  fluorescence  spectrometry  (AFS)  (1-10).  In 
general,  these  studies  have  revealed  that  laser-excited  AFS  can  provide  higher 
sensitivity,  reduced  interference  from  spectral  background  and  radiation 
scattering,  longer  linear  working  curves,  greater  freedom  from  quenching 
effects  and  source  fluctuations,  and  enhanced  multielement  capability  over 
AFS  excited  by  conventional  sources.  Several  of  these  advantages  derive 
directly  from  the  ability  of  high-power  tunable  lasers  to  saturate  an  atomic 
transition  (6,  11-13). 

bnergy  level  saturation  occurs  in  an  atomic  system  when  radiationally 
induced  excitation  becomes  much  more  rapid  than  all  deactivation  processes. 
Under  such  conditions,  the  state  from  which  the  excitation  transitions  ori- 
ginate becomes  somewhat  depopulated  while  the  excited  state  population  in- 
creases. Because  the  excitation  rate  producing  this  depopulation  depends 
upon  the  photon  flux  incident  on  the  absorbers,  the  upper  state  population 
can  approach  that  of  the  lower  state,  provided  that  the  source  irradiance  is 
sufficiently  high.  When  this  happens,  the  atomic  system  approaches  satura- 
tion and  is  said  to  be  bleached. 

If  an  atomic  system  is  bleached  or  saturated,  its  absorption  is  non- 
linear and,  in  fact,  approaches  zero,  since  the  atoms  cannot  absorb  more 
energy  than  is  necessary  to  produce  equivalent  upper  apd  lower  energy  level 
populations.  Furthermore,  as  long  as  the  incident  source  radiation  has  a 
sufficient  power  density,  saturation  will  be  maintained  despite  the  existence 
of  or  changes  in  radiationless  depopulation  of  the  excited  state.  Consequent- 
ly, the  detected  spontaneous  fluorescence  from  a saturated  atomic  system  will 
be  constant,  since  spontaneous  radiational  deactivation  rates  depend  only  on 
the  excited  state  population.  Thus,  atomic  fluorescence  measured  under  con- 
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ditions  of  high  source  irradiance  is  immune  to  source  power  variations,  as  long  as  the 
power  is  sufficient  to  maintain  saturation.  Also,  because  the  saturated  ex- 
cited state  population  is  governed  entirely  by  radiationally  induced  processes, 
quenching  of  individual  excited  atoms  by  environmental  (flame)  species  does 
not  significantly  affect  the  average  number  of  excited  atoms,  and  the  detected 
fluorescence  is  similarly  unaffected.  Working  curves  of  sample  concentration 
vs.  fluorescence  exhibit  extended  linearity  under  conditions  of  energy-level 
saturation,  due  to  the  absence  of  inner-filter  effects  (5-7,  10).  These  ef- 
fects become  unimportant  because  of  the  transparency  of  a "bleached"  atomic 
system  to  radiation  at  the  exciting  wavelength. 

Because  of  the  obvious  importance  of  energy-level  saturation  to  laser- 
excited  AFS,  it  would  be  useful  to  know  the  minimum  incident  radiant  power 
density  which  would  be  necessary  to  achieve  saturation  of  the  energy  levels 
in  an  atomic  system.  Such  knowledge  would  enable  one  to  take  full  advantage 
of  the  practical  characteristics  produced  by  the  saturation  effect  and  would 
allow  a prediction  of  the  feasibility  of  producing  saturation  by  specially 
designed  conventional  (i.e.  non-laser)  sources. 

Previous  research  work  in  saturation  of  atomic  energy  levels  achieved 
at  high  spectral  irradiance  has  predicted  little,  if  any,  effect  of  the 
quenching  characteristics  of  the  atom  environment  on  the  saturation  power 
density  (5,  6,  10-13).  Also,  all  the  theoretical  models  proposed  to  explain 
the  saturation  effect  have  assumed  that  steady-state  conditions  prevail  in 
the  radiationally  induced  excitation  process,  an  assumption  which  is  only 
valid  if  the  duration  of  the  excitation  pulse  is  much  longer  than  the  flu- 
orescence lifetime  of  the  transition  being  studied.  Because  typical  values 
of  atomic  fluorescence  lifetimes  are  1-20  ns,  such  steady  state  conditions 
do  not  prevail  in  atomic  systems  excited  by  a short-pulsed,  nitrogen-pumped 
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tunable  dye  laser  (typical  pulse  width:  5-10  ns),  although  they  are  approached 
when  longer  duration  sources  are  used.  One  such  source  is  the  flash- lamp- 
pumped  tunable  dye  laser  (typical  pulse  lengths  1-10  us) . 

In  the  present  study,  experimental  and  predicted  values  of  saturation 
power  density  were  determined  and  their  dependence  on  atom-environment  and 
on  the  spectral  and  spatial  characteristics  of  the  irradiance  source  was  as- 
certained. The  study  was  comprised  of  and  will  be  reported  in  two  parts. 

In  the  part  covered  here,  a theoretical  model  will  be  developed  for  saturation 
of  atomic  levels  by  a radiation  pulse  of  variable  duration.  The  model  includes 
provision  for  radiationless  (collisional)  deactivation  of  excited  species  and 
is  applicable  to  complex  radiational  decay  schemes  as  well  as  to  simpler  res- 
onance fluorescence  behavior.  To  illustrate  the  utility  of  the  model,  values 
of  the  radiant  spectral  power  density  necessary  to  saturate  specific  atomic 
transitions  will  be  calculated  and  presented,  and  the  dependence  of  saturation 
spectral  power  den-ity  on  excitation  pulse  width  will  be  illustrated. 

In  the  second  part  of  the  study,  to  be  published  later,  experimental 
values  of  saturation  spectral  power  density  will  be  reported  for  a number  of 
atomic  transitions,  using  a nitrogen-pumped  dye  laser  as  the  excitation 
source  (pulse  width  = 6.5  ns).  Conditions  necessary  for  obtaining  accurate 
values  of  saturation  spectral  power  density  will  be  described  and  the  influ- 
ence of  experimental  parameters  such  as  atom  environment,  atom  concentration, 
and  type  of  transition  will  be  illustrated.  Finally,  experimental  values  of 
saturation  spectral  power  density  so  obtained  will  be  compared  with  those 
predicted  by  the  model  developed  herein;  the  resulting  agreement,  considerably 
improved  over  that  obtainable  using  earlier  theoretical  models,  enables  useful 
conclusions  to  be  drawn  concerning  the  feasibility  and  practicality  of  employ- 
ing saturated  atomic  fluorescence  spectrometry  for  routine  analysis. 
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In  this  section,  a theoretical  model  is  proposed  to  explain  the  depen- 
dence of  the  non-steady- state  saturation  of  atomic  energy  levels  on  the  rela- 
tive width  of  the  excitation  pulse  and  on  the  fluorescence  lifetime  of  the 
observed  atomic  transition.  Also  considered  theoretically  are  the  influence 
of  quenching  species  in  the  atom's  environment  on  the  saturation  power  den- 
sity and  the  importance  of  non-resonance  transitions  to  the  saturation  process. 

Relationship  Between  Upper-State  Population  and  Fluorescence  Yield. 

In  general,  for  a spontaneous  atomic  fluorescence  transition  from  an 
upper  energy  level  2 of  energy  E2  to  a lower  energy  level  1 of  energy  Ei,  the 
magnitude  of  the  detected  fluorescence  radiance  (Bp)  at  a given  time  (t)  is 
given  by 

BF  (t)  = C A2 1 n2(t)  (E2  - EO  (1) 

where  C is  a proportionality  constant  that  contains  geometrical  parameters 
such  as  angle  of  collection  of  the  fluorescent  radiation,  the  length  of  the 
absorbing  atomic  medium,  etc.;  A2 1 is  the  Einstein  coefficient  for  the  spon- 
taneous transition  2 -*•  1;  n2(t)  is  the  atom  population  in  the  excited  energy 

level  2 at  time  t;  E2  - Ei  = hv2j  is  the  energy  difference  between  the 
levels  involved  in  the  transition  2 -*■  1,  where  h is  the  Planck  constant  and 
v2 1 the  spectral  frequency  of  the  2 -+  1 transition. 

From  equation  1,  the  detected  fluorescent  radiation  obtainable  from  any 
number  of  excited  atoms  is  dependent  only  on  the  optical  design  of  the  radia- 
tion-collecting system  being  used,  on  fundamental  constants  which  are  charac- 
teristic of  the  particular  transition  being  viewed  and  on  the  number  of  atoms 
that  exist  in  the  excited  state  from  which  the  observed  fluorescent  transition 
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originates.  Of  these  parameters,  only  the  excited  state  population  is  of 
significance,  the  others  usually  being  fixed  for  any  given  experimental 
arrangement . 

Because  the  detected  fluorescence  is  directly  proportional  to  the  time- 
dependent  excited-state  population  [n2(t)],  it  becomes  possible  to  study  the 
practical  significance  of  energy-level  saturation  in  atomic  fluorescence 
spectrometry  merely  by  examining  the  effects  of  saturation  on  the  population 
of  specific  upper  energy  levels  involved  in  the  desired  fluorescence  process. 
For  this  reason,  all  subsequent  treatment  of  saturation  processes  in  this 
chapter  will  be  oriented  toward  determining  particular  excited  state  popula- 
tions, with  the  understanding  that  such  populations  can  be  directly  related 
to  fluorescence  yields. 

Since  resonance  and  direct- line  fluorescence  (14)  are  the  most  frequently 
used  At-  transitions  in  chemical  analysis,  we  will  separately  consider  a two- 
level  and  a three- level  atomic  model,  each  yielding  one  of  these  two  kinds  of 
fluorescence  transitions.  These  models  will  then  be  used  to  explain  the  time 
behavior  of  the  population  of  the  excited  energy  level  involved  in  fluorescence 
as  a function  of  the  length  of  the  radiation  pulse  causing  excitation. 

A diagram  of  the  two  energy-level  (resonance)  atomic  model  and  the  math- 
ematical function  to  describe  the  time  behavior  of  the  exciting  radiation 
pulse  are  shown  in  Fig.  1.  Fig.  1A  illustrates  a resonance  transition  with 
upper  (2)  and  lower  (1)  atomic  levels  of  energies  E2  and  Ei  respectively, 
with  the  relevant  activation  and  deactivation  paths  of  the  transition  indicated 
by  arrows.  In  this  case,  the  atomic  system  is  assumed  to  be  irradiated  by  a 
beam  of  photons  having  a spectral  power  density  p(t)  which  is  constant  over 
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the  wavelength  (frequency)  region  of  interest  and  which  has  a time  variation 
described  by  a square  pulse  of  definable  length  (Fig,  IB). 

In  tli i s model,  a group  of  atoms  in  the  lower  energy  level  1,  when  exposed 
to  the  beam  of  radiation,  moves  to  the  upper  energy  level  2 by  absorption  of 
a quantum  each  of  the  incident  energy,  and  at  a rate  proportional  to  Bi2P(t)- 
B 1 2 is  the  Einstein  coefficient  for  the  stimulated  absorption  transition  1 2 . 

Once  in  the  excited  level,  the  atom  can  undergo  deactivation  to  the  lower  energy 
level  either  by  spontaneous  emission,  by  collision  with  other  species  or  through 
emission  induced  by  the  incident  radiation.  The  rate  constants  of  these  deac- 
tivation processes  are  (respectively): 

A2i(sec_1);  the  Einstein  coefficient  for  spontaneous  emission,  which 
indicates  the  probability  that  an  excited  atom  will  spontaneously  emit, 
in  a random  direction,  a quantum  of  energy  equal  to  that  of  the  ab- 
sorbed radiation  (E2  - Ei); 

Z2i(sec"1);  the  rate  of  collisional  deactivation,  or  the  probability 
that  an  excited  atom  will  pass  to  the  lower  energy  level  1 by  collision 
with  other  species  without  emission  of  radiation; 

b2  lp (t)  (sec-1);  the  probability  that  the  excited  atom,  exposed  to  the 
excitation  beam  of  power  density,  p(t),  will  return  to  the  lower  energy 
level  1 by  emission  of  a quantum  in  the  same  direction  as  the  in  ident 
radiation.  B2i  is  the  Einstein  coefficient  for  the  radiationally  in- 
duced downward  transition  2 -*■  1. 

This  latter  kind  of  deactivation  process  is  not  detected  in  conventional 
AFS,  even  though  it  is  radiational,  because  of  its  directional  character. 

B12  and  B2 1 are  related  to  A2  1 and  each  other  by 


2-^b21 
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and 

B21  *2  = B12  8i  (3) 

where  h is  the  Planck  constant  (6.624  * 10-34  J sec),  c the  speed  of  light 
(3  x lo10  cm/sec),  X21  is  the  wavelength  of  the  2 -*■  1 transition  and  gi,  g2 
are  the  degeneracies  (statistical  weights)  of  energy  levels  1 and  2 respectively. 

llius,  considering  only  the  activation  and  deactivation  processes  mentioned 
above,  and  neglecting  temperature  and  pressure  broadening  of  the  atomic  spectral 
lines,  the  rate  of  population  of  the  excited  energy  level  2 of  this  model  is 
given  by: 

~4-=  B12  p(t)  nx(t)  - (Bal  p(t)  + Aai  + Z2i)  n2(t)  (4) 

where  nj(t)  and  n2(t)  are  the  time-dependent  number  of  atoms  in  energy  levels 
1 and  2,  respectively. 


Assuming  that  the  total  number  of  atoms  (N)  is  distributed  among  those 
two  levels, 


N = t)  + t) 

(5) 

By  substitution  of  n 

j(t)  from  Eq.  (5)  into  Eq.  (4),  the  rate 

of  population 

excited  level  can  be 

found  in  terms  of  n2(t): 

NBi2P(t)  ■ (a2  + baP(t)} 

(6) 

where 

Agi  + Z21 

(7) 

and 

»2=  B12  + B21 

(8) 

It  is  noteworthy  that  the  constant  a2  is  the  reciprocal  of  the  fluorescence  life- 
time (x')  of  the  observed  transition. 

In  this  analysis,  the  time  behavior  of  the  laser  pulse  [p(t)]  is  assumed 
to  be  a boxcar  function  of  width  t0  and  height  p0  (cf.  Fig.  IB),  mathematically 

described  by 


*<«>  - { So 

0 


0 6 t S t0 


t > t 
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Remembering  the  shape  of  this  excitation  pulse,  let  us  first  consider  the  case 
in  which  a steady-state  population  of  the  excited  level  is  achieved  during  the 
excitation  pulse;  this  situation  prevails  when  the  excitation  pulse  length 
(to)  is  large  compared  to  the  fluorescence  lifetime  of  the  transition  (t'). 


Case  I:  Steady  State  Irradiation 


Under  steady-state  irradiation  conditions,  the  change  in  population  of  the 


excited  state  is  negligible  so  that 


dn2(t)  ^ n 

TZ 'U  U 


If  the  exciting  radiation  then  has  a constant  power  density  of  p q (i.e., 
t0  ► “) , the  steady-state  population  of  energy  level  2 [N§s],  obtained  by 


solving  Eq.  (6)  is  given  by 


= OT12  C 


12  V a2  + b2po 


where 


Ps  b . 


Eq.  (10)  implies  that  under  steady-state  conditions,  the  population  of 
the  excited  level  increases  in  nonlinear  proportion  with  the  peak  power  den- 
sity of  the  excitation  pulse  (po),  toward  a saturation  level,  [(N2)sat],  which 
can  be  found  from  Eq.  (10).  In  Eq.  (10),  as  po  gets  larger  and  larger,  ps 
becomes  less  and  less  significant.  Ultimately,  the  denominator  in  the  paren- 
thetical expression  of  Eq.  (10)  becomes  approximately  equal  to  p0  and  Eq.  (10) 
simplifies  to 


2 sat  b- 
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Thus,  under  steady-state  conditions,  the  saturation  population  of  an 
atomic  system  depends  only  on  the  fundamental  parameters  B12  and  b2  = ti\2  * Bj-  j 
[Eq . (8) j . From  Eq.  (3),  the  values  of  B12  and  B21  are  seen  to  be  related 
by  their  statistical  weights.  Using  this  information  and  Eq.  (12),  the  maxi- 
mum fraction  of  excited  atoms 


N vj.  ' ^ (13) 

is  only  dependent  on  the  degeneracies  of  the  levels  involved. 

In  practical  work,  it  has  been  found  useful  to  define  a saturation  power 
density  (7)  which  occurs  when  po  = pg.  In  this  situation,  Eq.  (10)  becomes 

NB,„ 

(NSS)  -r-2£(i) 

2 3 b2 

or 


(Os  - * (Vsat 


In  other  words,  when  the  incident  power  density  po  reaches  the  value  ps,  the 
excited-state  population  will  be  50  percent  of  that  achievable  under  complete 
saturation . 

Case  II.  Non-Steady  State  (Transient)  Irradiation 

When  the  duration  (to)  of  the  radiation  pulse  producing  excitation  is 
comparable  to  or  shorter  than  the  fluorescence  lifetime  of  the  upper  state 
of  the  transition  being  observed,  the  steady-state  assumption  made  in  the  pre- 
vious section  does  not  hold,  and  the  change  in  population  of  the  excited  level 
[dn2(t)/dt]  cannot  be  neglected.  In  this  non-steady-state  situation,  the  popu- 
lations of  the  excited  level  both  during  (0  .1  t i.  to)  and  after  (t  > to)  ter- 
mination of  the  excitation  pulse  can  be  obtained  by  solving  Eq.  (6),  and  are 
given  by: 
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A plot  of  the  above  expressions  reveals  the  time  behavior  of  the  excited 
state  population  expected  when  exciting  radiation  pulses  of  any  length  are 


employed.  For  convenience,  Fqs.  (14)  and  (15)  will  be  plotted  as  the  fraction 


of  excited  atoms  — — , rather  than  directly  as  upper-level  populations.  As 


an  example.  Fig.  2 shows  the  changing  fractional  upper-level  population  ex- 
pected under  conditions  prevailing  during  excitation  by  a nitrogen-pumped 
tunable  dye  laser.  To  construct  this  plot,  an  excitation  pulse  of  5 ns  dura- 
tion was  assumed;  this  value  is  typical  of  pulse  lengths  in  nitrogen-pumped 
tunable  dye  lasers.  For  the  > b2P0  and  B l 2 parameters,  values  typical  of 
many  atomic  fluorescence  transitions  were  used,  specifically; 


a =■  10' ' Sec’i 


b $ a 

Or.  O 

t > C. . 

II  --  b J'd.  (g  - g in  F.qs.  (3)  and  (8)) 

In  Fig.  2,  the  fraction  of  excited  atoms  rises  as  long  as  the  exciting 
radiation  is  applied;  the  rise  time  of  the  population  increase  is  calculated 
to  be  0.5  ns,  which  corresponds  to  the  quantity  (a2  + b2Po)-1>  as  one  would 
expect  from  Lq.  (14).  After  termination  of  the  excitation  pulse,  the  popula- 
tion of  the  excited  level  decays  with  a time  constant  of  1 ns,  which  corresponds 
to  the  effective  fluorescence  lifetime  of  the  atomic  transition  in  a flame 
environment  [ x ’ = aj1 , cf.  Kq.  (15)].  In  this  special  case,  the  half-width 
of  the  curve  describing  the  fraction  of  excited  atoms  is  equal  to  5 ns,  the 


r 
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same  as  the  length  of  the  excitation  pulse. 

As  mentioned  earlier,  the  variation  in  excited-state  population  illu- 
strated in  Fig.  2 can  be  experimentally  followed  by  monitoring  the  spontaneous 
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fluorescence  emitted  during  depopulation  of  the  excited  state.  Because  the 
magnitude  of  spontaneous  fluorescence  is  directly  proportional  to  the  excited- 
state  population  [cf.  Eq.  (l)],a  plot  of  measured  fluorescence  vs.  time  should 
generate  a curve  such  as  Fig.  2 directly.  However,  because  of  limitations 
in  instrumentation,  it  is  often  more  convenient  to  measure  the  average  fluo- 
rescence which  is  produced  over  a specified  time  interval.  In  the  specific 
but  frequently  encountered  example  of  gated  detector  (boxcar  integrator)  being 
employed,  the  recorded  signal  will  be  the  time  average  of  the  fluorescence 
radiation,  detected  over  the  sampling  gatewidth  of  the  gated  detector.  For 
convenience,  let  us  assume  a sampling  gatewidth  equal  to  the  laser  pulse  dura- 
tion (to).  The  average  population  of  the  excited  energy  level  (N2)  during 
the  excitation  pulse  (which  determines  the  magnitude  of  the  detected  fluores- 
cence signal)  can  then  be  calculated: 

__  _1 

N = n ( t)  = t 

2 2 o 

By  substitution  of  Eq.  (14)  into  Eq.  (16)  and  integration,  one  obtains 

"z  - miz  (aa  Asp0  ) n - P)  (17) 

where 

D = 


As  in  the  case  of  steady-state  population  of  the  excited  state,  a value  for 
the  radiant  power  density  can  be  defined  (p^^)  which  is  necessary  to  produce 
an  upper  state  population  just  half  of  that  required  for  saturation.  Thus, 


_ e~(a2  + b2P<j)  tp) 
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by  solving  Eq.  (17)  for  po  = Pg  ^ and  doing  a series  expansion  of  the  expo- 
nential term  in  the  expression  for  I),  one  obtains  the  following  expression 


for 


NSS 

3s  - 


(b2tQ) 


-l 


(19) 


wh i oh  implies  that  under  non-steady-state  population  of  the  excited  level, 
the  saturation  spectral  power  density  depends  only  on  the  Einstein  coeffi- 
cients for  induced  transitions  ( B 1 2 » B2 l ) and  on  the  duration  of  the  excita- 
tion pulse  but  is  independent  of  the  rate  of  quenching  collisions.  A com- 
parison of  Eqs.  (10)  and  (17)  reveals  that  under  non-steady-state  conditions, 
the  rate-controlled  average  population  of  the  excited  level  (N2)  differs  from 
the  steady-state  population  (Nfs)  by  an  amount  given  by  the  D term  (cf.  Eq. 
18). 

According  to  Eq.  (18),  the  departure  (D)  of  N2  from  is  determined, 
at  any  given  value  of  p0,  by  the  relative  magnitudes  of  the  excitation  pulse 
length  (to)  and  the  excitation  time  constant  given  by  (02  + b2Po)_1-  For 
values  of  to  much  larger  than  (a  + bpo)-1,  f>  becomes  negligible,  and  N2  be- 
comes equal  to  n|s,  an  altogether  expected  result.  In  contrast,  at  values 
of  to  comparable  to  or  smaller  than  (a2  + b2Po)_1,  the  population  of  the  ex- 
cited level  (N2)  rises  with  increasing  po  toward  the  saturation  level  at  a 
rate  determined  primarily  by  the  term  D and,  therefore,  by  the  to,  a and  b 
parameters. 

The  rate-controlled  population  process  described  above  is  illustrated 
in  Fig.  3,  for  excitation  pulses  of  several  different  lengths.  In  Fig.  3, 
the  fraction  of  excited  atoms  (N/N)  has  been  calculated  from  Eq.  (17)  and 
is  plotted  as  a function  of  the  excitation  power  density  pq  (watts/cm2  Hz), 
for  excitation  pulses  of  l ns,  10  ns  and  1 ps  duration.  These  particular 
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plots  correspond  to  the  behavior  predicted  for  the  So  + Pi  transition  of  Ba 
(553.5  mi) . For  this  transition,  a fluorescence  lifetime  (t)  of  8.5  ns  has 
been  measured  (15)  in  the  absence  of  foreign  quenching  gases  (i.e.,  when 
Z2 i ' 0).  For  the  plots  in  Fig.  3,  the  a2  and  b2  parameters  were  calculated 
by  assuming  that  the  rate  of  collisional  deactivation  would, in  a flame  en- 
vironment, be  approximately  ten  times  faster  than  the  rate  of  deactivation  by 
spontaneous  emission  (i.e.,  Z2 1 % 10  A21)  and  by  using  the  relationship  be- 
tween t)ie  A2iand  B2 1 coefficients  given  in  Eqs.  (2)  and  (3). 

From  Fig.  3,  the  population  of  the  excited  level  rises  toward  the  satur- 
ation level  faster  when  a 1 nsec  pulse  is  used  (curve  A)  than  when  a 1 nsec 
excitation  pulse  is  employed  (curve  C) . These  results  are  entirely  expected; 
when  the  excitation  pulse  (to)  is  much  longer  (1  usee)  than  the  fluorescence 
lifetime  of  the  transition  and  therefore  also  much  longer  than  (a2  + b2Po)_1j 
steady-state-population  of  the  excited  level  should  prevail.  In  contrast,  a 
1 ns  duration  excitation  pulse  requires  non-steady- state  (kinetically  controlled) 
population  of  the  excited  level,  because  the  length  of  the  pulse  is  smaller 
than  the  fluorescence  lifetime  of  the  transition  and  comparable  to  (a2  + b2Po)_: 
for  low  values  of  Po.  When  a 10  ns  excitation  pulse  is  assumed  (curve  B) , 
the  duration  of  the  pulse  is  comparable  to  the  fluorescence  lifetime,  but 
longer  than  (a2  + b2Po)-1.  Therefore,  the  fraction  of  excited  atoms  varies 
with  exciting  power  density  just  as  in  curve  A,  but  with  a lower  magnitude; 
this  behavior  arises  because  the  D factor  (cf.  Eq.  (18))  is  not  negligible 
in  this  case. 

A comparison  of  curves  A and  C in  Fig.  3 reveals  that  the  extent  of 
departure  of  eacli  curve  from  the  behavior  expected  for  steady-state  popula- 
tion decreases  as  the  spectral  power  density  po  increases.  In  fact,  at  very 
high  values  of  p0,  the  fraction  of  excited  atoms  should  become  the  same  in  all 
cases. 


1 


14 


In  order  to  calculate  the  saturation  power  density  corresponding  to 
each  of  the  curves  in  Fig.  3,  we  must  know  the  maximum  fraction  of  excited 
atoms  which  would  exist  during  complete  saturation  of  the  atomic  system.  In 
the  case  of  the  So  Pi  transition  of  Ba  under  consideration,  a maximum  frac- 
tion of  excited  atoms  equal  to  0.75  was  calculated  from  Eq.  (13),  with  gi  = 1 
and  g;.  = 3.  From  Fig.  3,  the  saturation  power  density  (i.e.,  the  spectral 
power  density  necessary  to  produce  a fraction  of  excited  atoms  equal  to  0.375) 
varies  with  the  duration  of  the  excitation  pulse. 

Table  I shows  calculated  values  of  saturation  power  densities  for  the 
excitation  pulses  considered  in  Fig.  3.  From  Table  I,  a 1 ns  exciting  pulse 
must  be  over  twice  the  power  density  of  a 1 ps  pulse  to  produce  an  upper 
level  population  equal  to  50  percent  of  that  prevailing  at  total  saturation. 
This  result  is  a consequence  of  the  kinetically  controlled  nature  of  the 
saturation  process  which  must  be  dominant  when  short  exciting  pulses  are  used. 

From  Eq.  (17),  the  dependence  of  the  excited-state  population  on  the 
excitation  power  density  is  determined  by  the  magnitude  of  the  a2  and  b2 
lumped  parameters  for  a given  excitation  pulse  duration.  Therefore,  the 
saturation  power  density  which  is  predicted  by  the  described  model  is  deter- 
mined by  both  fundamental  atomic  parameters  and  environmental  factors,  such 
as  the  fluorescence  lifetime  of  the  transition  and  the  rate  of  collisional 
deactivation  produced  by  the  atom  environment  (cf.  Eqs.  (7)  and  (8)).  Un- 
fortunately, accurate  values  for  collisional  deactivation  rates  in  a flame 
environment  are  not  often  available;  furthermore,  the  fluorescence  life- 
times listed  for  most  atomic  transitions  contain  significant  uncertainties. 

To  better  appreciate  the  utility  of  this  theoretical  approach,  let  us  then 
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consider  the  effect  of  the  accuracy  of  these  parameters  on  the  predicted 
saturation  curves  and  on  the  calculated  saturation  power  densities.  In  the 
following  considerations,  the  same  transition  of  Ba  which  was  described 
earlier  in  this  section  will  be  employed. 

First,  let  us  assess  the  importance  of  uncertainties  in  the  rates  of 
collisional  relaxation  (Z2i).  In  a flame  or  other  atomizer  environment,  the 
fluorescence  lifetime  of  any  given  atomic  transition  will  be  shortened  by 
collisions  between  excited  atoms  and  other  flame  species.  Such  behavior 
was  verified  experimentally  in  this  laboratory  by  measurement  of  fluorescence 
lifetimes  of  atomic  transitions  in  the  flame  environment  and  comparing  those 
values  to  lifetimes  of  isolated  atoms. 


'Hie  effect  of  the  rate  of  collisional  deactivation  upon  theoretical 
saturation  curves  is  shown  in  Figs.  4 and  5 for  excitation  pulse  lengths  of 
1 gs  and  6.5  ns  respectively.  In  constructing  Fig.  5,  we  have  assumed  a 
6.5  ns  pulse  length  because  it  corresponds  to  the  excitation  pulse  of  a 
n itrogen- laser-pumped  tunable  dye  laser  which  was  employed  in  the  experimental 
measurements  of  saturation  power  densities  to  be  described  in  a later  commu- 
nication. In  both  Figs.  4 and  5,  curve  A corresponds  to  the  saturation  curve 
where  the  atom  environment  is  free  of  foreign  quenching  species  (Z2i  ^ 0), 
while  curves  B and  C correspond  to  atom  environments  where  the  rate  of  col- 
lisional deactivation  is  assumed  to  be  5 and  10  times  faster  than  the  rate 


of  spontaneous  emission  (Z2i  = 5A2i,  IOA21),  respectively. 

When  a 1 ps  excitation  pulse  is  used  to  excite  an  atomic  medium  rela- 
tively free  of  quenching  species  (Fig.  4A) , saturation  of  the  excited  level 

4 * 

is  achieved  at  relatively  low  values  of  incident  spectral  power  density. 
However,  in  the  presence  of  quenching  species  (Figs.  4B  and  4C) , the  fraction 
of  excited  atoms  rises  toward  the  saturation  level  at  a rate  which  is  in- 
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versely  proportional  to  the  rate  of  collisional  deactivation. 

If  a 6.5  ns  excitation  pulse  is  employed  (cf.  Fig.  5),  the  saturation 
curve  is  affected  much  in  the  same  way  by  collisional  de-excitation  as  when 
the  longer  (1  us)  pulse  is  used.  However,  because  the  excited  state  popu- 
lation is  determined  mainly  by  radiationally  induced  transitions  when  the 
o.5  ns  pulse  is  applied,  saturation  is  achieved  only  at  higher  spectral 
power  densities  than  for  the  1 us  pulse,  where  a steady-state  population 
of  the  excited  level  prevails.  A comparison  of  curves  A,  B,  and  C in  either 
Figs.  4 or  5 reveals  that  the  atomic  system  approaches  saturation  as  the 
spectral  power  density  increases  regardless  of  the  occurrence  of  quenching 
collisions.  In  fact,  at  spectral  power  densities  much  higher  than  the  sat- 
uration power  density  (ps),  the  fraction  of  excited  atoms  becomes  essentially 
independent  of  the  rate  of  collisional  deactivation.  This  behavior  is  ex- 
pected because  at  high  power  densities,  radiationally  induced  transitions 
(b2p 0 ) become  dominant  in  the  population  and  depopulation  of  the  excited 
energy  level. 

The  calculated  saturation  spectral  power  densities  corresponding  to 
the  different  curves  in  Figs.  4 and  5 are  given  in  Table  II.  In  this  table, 
t is  the  fractional  (expressed  as  a percentage)  difference  between  the  sat- 
uration power  densities  which  were  calculated  for  the  6.5  ns  and  1 usee  ex- 
citation pulse  lengths  (columns  2 and  3 of  Table  II).  From  a practical 
standpoint,  c expresses  the  relative  error  which  would  exist  in  the  calcu- 
lated value  of  saturation  power  density  if  steady-state  conditions  had  been 
assumed,  but  a 0.5  ns  excitation  pulse  were  employed. 

Several  interesting  and  important  observations  can  be  made  from  Table 
II.  When  a 1 us  excitation  pulse  is  assumed,  there  is  a one  order  of  mag- 
nitude difference  between  the  values  of  calculated  for  an  atomic  medium 
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relatively  free  of  quenching  species  (Z21/A21  £ 0)  and  that  pertaining  to 
a medium  where  the  rate  of  collisions  is  one  order  of  magnitude  greater 
than  the  rate  of  spontaneous  emission  (Z21/A21  = 10).  In  comparison,  the 
dependence  of  P s on  quenching  is  smaller  when  a 6.5  ns  excitation  pulse  is 
assumed,  although  the  values  of  Ps  are  higher  than  those  required  for  the 
1 us  pulse.  Also,  Table  II  indicates  that  the  error  produced  by  assuming 
a steady-state  population  of  the  excited  level,  in  the  case  of  a 6.5  ns 
pulse,  decreases  as  the  quenching  factor  (Z21/A21)  increases.  This  behavior 
arises  because,  as  collisional  deactivation  (Z21)  increases,  the  effective 
fluorescence  decay  rate  of  the  atomic  transition  (a2)  becomes  comparable  to 
or  smaller  than  the  duration  of  the  excitation  pulse.  Consequently,  the 
6.5  ns  pulse  begins  to  appear  more  and  more  like  a continuous  (in  time) 
radiation  level,  so  that  steady-state  irradiance  conditions  are  approached. 

Another  parameter  which  might  cause  significant  error  in  the  calculated 
values  of  saturation  power  density  is  the  fluorescence  lifetime  of  the  tran- 
sition being  considered.  The  uncertainties  in  published  lifetimes  are  sub- 
stantial, and  the  value  finally  used  in  the  calculation  directly  determines 
the  rate  of  depopulation  of  the  excited  energy  level,  thereby  strongly  af- 
fecting the  incident  power  density  required  to  produce  saturation.  A plot  of 
saturation  curves  obtained  from  Eq.  (17),  assuming  values  of  4,  8 and  16  ns 

for  the  fluorescence  lifetime  of  the  Ba  transition  considered  and  a 6.5  ns 


pulse  duration  are  shown  in  curves  A,  B and  C of  Fig.  6.  As  expected,  the 
fraction  of  excited  atoms  is  greatly  dependent  upon  the  assumed  value  of  the 

i 

fluorescence  lifetime.  Predicted  values  of  saturation  power  density  corres- 
ponding to  the  curves  in  Fig.  6 are  shown  in  Table  III.  In  this  table,  the 
values  of  ps  increase  as  t becomes  larger  than  the  duration  of  the  excitation 
pulse.  This  table  also  shows  the  error  (e)  in  the  predicted  saturation 
power  density  which  would  be  caused  by  assuming  values  of  the  fluorescence 
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lifetime  that  are  50  percent  smaller  or  higher  than  the  "real"  fluorescence 
lifetime  (8  ns)  of  the  studied  transition. 

The  foregoing  considerations  reveal  that  an  accurate  knowledge  of  the 
fluorescence  lifetime  of  the  transition  and  of  the  rate  of  collisions  is 
necessary  in  order  to  predict  accurate  saturation  power  densities  of  an 
atomic  system  irradiated  with  a pulse  of  a given  duration.  Using  the  two- 
energy- level  model  described  earlier,  predicted  values  of  the  saturation 
power  density  for  several  elements  are  shown  in  Table  IV.  In  the  table  is 
also  listed  the  wavelength  of  the  particular  transition  for  each  element 
being  considered,  and  the  values  of  the  fluorescence  lifetime  available  in 
the  literature.  In  calculating  these  saturation  power  densities,  we  have 
assumed  the  rate  of  collisions  to  be  ten  times  that  of  spontaneous  emission 
(Z21  % 10A21)  and  an  excitation  pulse  of  6.5  ns  duration. 

The  effect  of  the  assumed  fluorescence  lifetime  on  the  saturation  power 
density  of  the  excited  level  becomes  particularly  strong  in  the  case  of  Cd 
which  lias  a relatively  long  fluorescence  lifetime  (35  ns).  Table  IV  also 
shows  that  a higher  spectral  power  density  is  required  to  saturate  atomic 
levels  in  the  ultraviolet  region  (e.g..  Mg)  compared  to  atomic  levels  in  the 
visible  spectral  region  (e.g.,  Ca,  In,  Sr,  etc.).  This  behavior  is  consis- 
tent with  the  dependence  of  the  Kinstein  A coefficient  for  spontaneous  emis- 
sion on  the  cube  of  the  frequency  of  the  transition  being  studied. 

It  is  both  worthwhile  and  instructive  to  consider  the  limits  of  the 


equations  describing  steady-state  and  non-steady-state  population  of  an  ex- 
cited level  [ bqs . (10)  and  (17)]  under  conditions  of  extremely  low  or  high 
spectral  irradiance.  These  special  cases  will  be  examined  separately. 
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Limits  Occurring  During  Steady-State  Excitation 

ior  spectral  power  densities  (pq)  much  less  than  the  saturation  power 
density  (pq  <<  Ps)>  the  population  of  the  excited  level  becomes,  from  Eq.  (10): 
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Eq.  (20)  predicts  that  under  steady-state  illumination  at  low  spectral  power 
densities,  the  population  of  an  excited  level  will  be  linearly  proportional 
to  the  effective  fluorescence  lifetime  of  the  transition  (t*  = a^1)  and  will 
also  increase  in  linear  proportion  with  the  power  density  (po)  of  the  excita- 
tion source.  This  is  just  the  behavior  observed  in  conventional -source  atomic 
fluorescence . 

At  excitation  source  densities  much  higher  than  the  saturation  power 
density  (ps)  the  population  of  the  excited  level  will  become  independent  of 
the  spectral  radiance  of  the  source  and  of  the  fluorescence  lifetime  of  the 
transition  as  given  by  Eqs.  (12)  and  (13).  This  behavior  has  been  observed 
by  other  workers  (6) . 


Limits  Occurring  During  Non-Steady-State  Excitation 

At  low  spectral  irradiances,  the  probability  of  induced  transitions 
(b2Po)  becomes  negligible  compared  to  the  probabilities  of  collisional  deac- 
tivation and  spontaneous  emission  (namely,  b2Po  <<c  ^2) . Therefore,  at  low 
incident  spectral  power  densities,  non-steady  state  population  of  an  excited 
level,  given  by  Eq.  (17),  becomes 
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From  Eq.  (21)  we  see  that  N2  (and  thus  any  detected  spontaneous  fluores- 


cence) increases  in  linear  proportion  with  pq,  but  its  magnitude  depends  on 
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the  relative  length  of  the  excitation  pulse  to  and  the  effective  fluores- 
cence lifetime  of  the  transition  t'  (t'  = aj1).  In  the  case  of  a source  of 

«•> 

low  spectral  radiance  and  pulse  duration  much  longer  than  x'  (conventional 
excitation  source),  the  population  of  the  excited  level  (Eq.  (21))  becomes 
equal  to  the  steady-state  population  at  low  spectral  irradiances  given  by 
Eq.  (20),  as  one  would  expect. 

At  luj^h  spectral  irradiances,  the  probability  of  radiational ly  induced 
transitions  (b^Po)  becomes  larger  than  the  probabilities  of  spontaneous  and 
collisional  deactivation  (a2).  Under  these  conditions,  Eq.  (17)  reduces  to 
the  following  relation 
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indicating  that  the  population  of  the  excited  level  becomes  independent  of 
the  effective  fluorescence  lifetime  x'  and  varies  in  nonlinear  proportion 
with  both  the  rate  of  radiational ly  induced  transitions  and  the  duration  of 

the  excitation  pulse. 

From  Eq.  (22),  for  an  excitation  pulse  of  longer  duration  than  (b2Po)_1 
or  for  a pulse  of  amplitude  (p0)  greater  than  ^tg)-1,  the  population  of 
the  excited  state  becomes: 
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which  is  the  same  as  the  maximum  population  of  the  excited  level  under  steady- 

0 

state  saturation  conditions  [Eqs.  (12)  and  (13)].  Eq.  (23)  implies  that  at 
high  power  densities  (p0  >>  ps)  or  under  steady-state  conditions  (to  >> 
(b2Po)-1),  the  atomic  medium  becomes  transparent  to  the  exciting  radiation 
and  the  population  of  the  excited  level  is  only  dependent  on  the  degeneracies 
of  the  energy  levels  involved  in  the  transition  and  on  the  atomic  population 
free  to  absorb  the  incident  radiation. 
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Because  of  the  high  available  peak  power  of  tunable  dye  lasers,  non- 
resonance transitions,  particularly  direct- lire  fluorescence  and  step-wise 
t iuorescence , have  become  potentially  useful  for  trace  analysis.  Therefore, 
it  is  important  to  consider  the  behavior  and  practical  applications  of  such 
transitions  under  conditions  of  near  saturation  of  the  relevent  atomic  en- 
ergy levels.  For  this  purpose,  we  will  consider  the  three-energy-level  model 
illustrated  in  Fig.  7A.  As  with  the  two-state  model,  it  will  be  assumed 
that  a number  (N)  of  three-state  atoms  is  illuminated  by  an  excitation  pulse 
resembling  a boxcar  function  of  duration  tg  and  peak  spectral  power  density 
P 0 . Such  a pulse  is  portrayed  in  Fig.  7B.  As  before,  the  exciting  radiation 
will  be  assumed  to  be  constant  in  magnitude  over  the  atomic  absorption  band- 
width. 

The  activation  and  deactivation  pathways  in  the  three-level  atom,  shown 
by  arrows  in  Fig.  7A,  correspond  to  the  following. 

Population  of  Level  2.  Excitation  of  an  atom  from  energy  level  1 to 
energy  level  2 can  occur  by  absorption  of  a quantum  of  incident  radiation. 

The  rate  of  this  process  is  given  by  Bj2  p(t),  where  B12  is  the  Einstein  co- 
efficient for  the  radiationally  induced  transition  1 -+■  2,  as  defined  earlier. 
Also,  an  atom  in  energy  level  3 can  be  excited  to  level  2 by  collision  with 
flame  species  or  through  other  thermal  processes.  The  rate  of  this  excita- 
tion event  is  given  by  Z32  and  becomes  important  only  in  cases  where  the 
energy  separation  between  energy  levels  2 and  3 is  comparable  to  or  smaller 
than  the  thermal  energy  of  the  flame  (i.e.,  E2  - E3  £ kT) . 

An  atom  originally  in  level  1 can  also  be  thermally  or  collisionally 
excited  to  level  2.  However,  the  energy  difference  between  these  states 
will  be  larger  than  that  between  levels  2 and  3,  so  that  such  collisional 
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excitation  is  less  probable.  In  the  present  case,  non-radiational  excita- 
tion from  level  1 to  level  3 will  be  neglected. 

Depopulation  of  Level  2.  Once  in  excited  level  2,  the  atom  can  relax 
through  the  following  paths: 

i)  by  radiational  deactivation  to  level  1,  involving  spontaneous  emission 

(in  a random  direction)  of  a quantum  of  energy  equal  to  that  in  the  inci- 
dent radiation.  This  process  is  called  resonance  fluorescence  and  its 
rate  or  probability  is  given  by  the  Einstein  coefficient  for  spontaneous 
emission  A2 1 . 

ii)  radiational  deactivation  to  level  3,  by  spontaneous  emission  (again  in 
a random  direction),  of  a quantum  of  energy,  (E?  - E3),  different  from 
that  of  the  incident  radiation.  The  rate  of  this  process,  termed  non- 
resonance direct-line  fluorescence,  is  given  by  the  Einstein  coefficient 
for  spontaneous  emission  A23. 

iii)  radiationless  deactivation  to  either  level  1 or  level  3 through  collisions 
with  foreign  species  in  the  flame.  The  rates  of  these  collisionally  in- 
duced transitions  are  given  by  Z21  and  Z23,  respectively, 

iv)  radiationally  induced  deactivation  to  level  1 with  emission  of  a quantum 
of  the  same  energy  and  in  the  same  direction  as  the  incident  radiation. 

The  rate  of  this  process  is  given  by  B21  p(t),  where  B21  is  the  Einstein 
coefficient  for  the  induced  transition  2 1. 

Population  and  Depopulation  of  Level  3.  Atoms  in  energy  level  3 will 
relax  to  energy  level  1 by  spontaneous  emission  at  a rate  given  by  A31  and 
by  collisions  with  flame  species  at  a rate  given  by  Z 3 1 . Excitation  of 
atoms  from  level  1 to  level  3 is  produced  by  collision  of  atoms  in  level  1 
with  flame  species  and  has  a rate  Z 3 3 ; collisional  excitation  from  level  3 
to  level  2 can  also  occur  and  will  have  a rate  Z32. 
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In  accordance  with  these  activation  and  deactivation  pathways,  the  rate 
of  population  of  energy  levels  2 and  3 can  be  shown  to  be: 


lt  ~ - n2^1p(t)  + a.i3  + a21j  + n3Z„ 


dn„( t) 
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and  nj,  nz  and  n3  are  the  time-dependent  atom  populations  of  energy  levels  1,  2 and  3, 
respectively.  In  this  and  the  following  treatment,  the  simplified  notation  n^  will  be 
used  rather  than  n^(t).  The  reader  should  recall  the  time-dependent  nature  of  these 
parameters,  however. 

Assuming  that  the  total  atom  population  (N)  is  entirely  distributed 
among  the  three  energy  levels  considered, 

N = ni  + n2  + n3  (29) 

In  the  case  of  direct-line  fluorescence,  the  signal  detected  is  the  fluores- 
cence radiation  produced  by  the  spontaneous  transition  from  level  2 to  level 
3,  Therefore,  the  pertinent  energy  level  in  determining  saturation  will  be 
upper  level  2.  By  contrast,  during  stepwise  fluorescence,  the  signal  detec- 
ted is  the  radiation  produced  by  the  spontaneous  transition  from  level  3 to 
level  1;  consequently,  the  population  of  energy  level  3 will  determine  the 
approach  to  saturation.  Further  discussions  of  saturation  of  this  three- 
level  atom  need  therefore  concern  only  the  populations  of  level  2 or  level 
3,  depending  on  whether  direct- line  or  stepwise  fluorescence  is  being  con- 
sidered . 
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The  probability  of  observing  either  direct-line  or  stepwise  fluores- 
cence, or  of  populating  any  given  energy  level  depicted  in  Fig.  7,  depends 
largely  on  the  rates  of  populating  or  depopulating  level  3.  To  see  why  this 
is  so,  let  us  consider  the  cases  of  direct-line  fluorescence  and  stepwise 

fluorescence  separately. 


Direct  1. 1 tie  Fluorescence 


Because  direct-line  fluorescence  ordinarily  involves  the  observation 
of  transitions  corresponding  to  frequencies  in  the  ultraviolet-visible  spec- 
tral region,  the  spacing  between  energy  levels  2 and  3 must  be  quite  large, 
and  will  generally  be  much  greater  than  the  thermal  energy  (kT)  available 
in  the  flame.  Consequently,  collisionally  induced  transitions  from  level  3 
to  level  2 will  generally  be  quite  slow  and  can  be  neglected  compared  to 
competing  processes  (i.e.,  Z32  0)  • This  assumption  will  be  valid  for 

some  atomic  systems  to  be  included  in  the  experimental  study  to  be  reported 
later.  For  other  systems  which  do  not  allow  this  assumption,  experimental 
correction  for  thermally  assisted  population  of  level  3 could  be  made  merely 
by  monitoring  emission  from  that  level  in  the  absence  of  exciting  radiation 
(i.e.,  when  p0  =0).  If  Z32  = 0,  the  rate  of  population  of  energy  levels  2 
and  3 [cf.  Eqs.  (24)  and  (25)]  becomes : 
dn2  (t) 

“ = niB1?p(t)  - n2(B2ip(t)  + a23  + a2i)  (30) 

and 


dn3( t) 
dt 


na  — n a 

2 23  1 13  3 31 
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To  study  the  saturation  of  upper  energy  level  2,  one  must  solve  Eq.  (30) 


i 


A 


as  a function  of  time  and  incident  power  density.  In  turn,  because  the  pop- 
ulation of  atoms  free  to  absorb  the  exciting  radiation  (namely  nj)  depends 
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partially  upon  the  rate  of  population  and  depopulation  of  energy  level  3, 
some  considerations  of  level  3 are  necessary  to  solve  Eq.  (30). 

Thermal  Equilibrium.  First,  consider  the  case  where  the  energy  differ- 
ence between  levels  3 and  1,  (E3  - Ej ) is  comparable  to  or  smaller  than  the 
thermal  energy  (kT)  available  from  the  flame  environment.  In  this  case,  an 
equilibrium  will  exist  between  the  population  of  levels  3 and  1 which  can 
be  described  by  the  relation: 


n Z « n„a^ 


Substitution  of  Eq.  (32)  into  Eq.  (29)  yields  the  population  of  level  1: 

(n  ' n2> 

n ^ ( 7 


where 


(1  + z) 


Z * Z1 3/  *31 


After  substitution  of  Eq.  (33)  into  Eq.  (30),  the  rate  of  population  of 


level  2 becomes: 
drt,(t)  Np(  t) B 


* "2  a21  + a; 


‘23  + [B21+  (1  + Z)J  } 
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where  the  lumped  parameters  are  defined  as: 


*3e  " “21  “23 


b = B2J  + 

3e 


1 (l  + Z) 


c3e  TTTzT 


In  these  expressions,  tiie  subscripts  3 and  e are  used  to  indicate  the  three- 
level  model  and  the  assumption  of  equilibrium  between  energy  levels  3 and  1. 
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A comparison  of  Eqs.  (6)  and  (35)  indicates  that  the  same  general  ex- 
pression for  the  rate  of  population  of  the  upper  level  is  obtained  for  a 
two-energy- level  atomic  system  and  a thrce-energy-level  system,  if  in  the 
latter  model  equilibrium  is  assumed  between  levels  3 and  1.  Also,  Eq.  (35) 
reduces  to  Eq.  (6)  for  the  case  Z % 0.  In  turn,  Z would  approach  zero  if 
energy  level  3 had  a negligible  population  compared  to  that  of  level  1 (cf. 
Eq.  34).  This  behavior  would  exist,  for  example,  in  the  case  of  a very  slow 
radiative  population  of  level  3 from  level  2,  coupled  with  a much  more 
rapid  radiative  depopulation  of  level  3.  In  such  a situation,  the  dominant 
pathway  for  energy  loss  from  level  2 would  be  the  transition  to  level  1, 
(i.e.,  A23  <<  A21).  In  such  a case,  one  need  consider  only  the  2 -*  1 de- 
excitation path  and  can  assume  a two-energy  level  system. 

Because  of  the  similarity  of  Eqs  (6)  and  (35),  the  same  analysis  which 
was  performed  for  the  two- level  system  on  the  population  of  the  upper  energy 
level  could  be  applied  to  the  present  case  of  a three  level  system,  with 
equilibrium  between  the  two  lowest  energy  levels  (3  and  1). 

Correspondingly,  a saturation  power  density  (P3es)  for  direct-line  flu- 
orescence can  be  defined  by  reference  to  Eq.  (11)  as 


P 


3ea 


(39) 
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By  substitution  of  Eq.  (40)  into  Eq.  (30),  one  obtains  the  following 
expression  for  the  rate  of  population  of  level  2: 
dn0( t)  NBt2p( t) 


dt 


— ~ (l — n z ) _ |a2i  + a23+  Beip(  t)  + B12o(  t) 


(l+ag3/ a31 ) 

(1  + z) 
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or 


dn  ( t)  r 1 

— » C38  P(t)  ' nM  [Ss  + 

dt 
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where : 


a3S  = + *23 


(1  + WSl) 


and 


b = B + B , . 

3*  21  12  (l  + Z) 

NB12 


(42) 


(43) 


C38  = (1  +“z7  (44) 

As  in  the  earlier  case,  the  population  of  the  excited  level  follows  the  same 
behavior  as  in  the  two-energy- level  system,  but  with  different  a,  b and  c 

coefficients  (cf.  Eq.  6). 


Stepwise  fluorescence 

In  atomic  systems  which  exhibit  stepwise  fluorescence,  the  energy  sep- 
aration between  excited  energy  levels  2 and  3 is  often  small  enough  that 
col  1 isional ly  induced  transitions  from  level  3 to  level  2 (Z32)  cannot  be 
neglected.  However,  the  rates  of  collisional  level- 1-to- level- 3 transitions 
become  negligible  because  of  the  large  energy  spacing  between  those  levels 
(corresponding  to  ultraviolet-visible  frequencies).  Under  these  conditions, 
Z13  ^ 0 and  the  rate  of  population  of  levels  2 and  3 are  (from  Eqs.  (24)  and 
(25) : 
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(45) 


and 


dn3( t) 


dt 


n2a23  - n'^a31  + Z32> 


(46) 


llecause  the  population  of  level  3 is  determined  solely  by  collisional 
and  non-radiationally  induced  transitions,  its  population  will  vary  as  a 
direct  function  of  that  in  level  2.  Likewise,  stepwise  fluorescence  detected 
as  a transition  from  level  3 to  level  1 will  follow  the  same  saturation  be- 
havior as  the  population  of  level  2. 

To  find  a simple  relationship  expressing  the  population  of  level  2,  it 
is  necessary  to  make  some  simplifying  assumptions.  One  such  assumption  is 
that  of  a steady-state  population  of  level  3.  If  the  rate  of  collisional 
deactivation  of  excited  atoms  from  level  2 to  level  3 is  comparable  to  the 
radiational  and  collisional  deactivations  of  level  3 to  level  1,  and  from 
level  3 to  level  2,  one  can  assume  a steady-state  population  of  level  3 (i.e, 
dn3(t)/dt  ^ 0).  In  this  case,  the  population  of  level  2 is: 


foaliQ  = NBiap(0  - ng(t){b  p(  t)  + a } 
dt  nr 


nr 
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a _ - a + a 
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(49) 


As  in  the  two  cases  treated  earlier  involving  direct-line  fluorescence,  the 
population  of  excited  level  2 follows  the  same  qualitative  dependence  on  the 
incident  power  density  as  did  the  upper  level  in  a resonance  transition. 
However,  in  the  present  case,  the  saturation  power  density  is  determined  by 
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the  rates  of  deactivation  from  level  2 to  level  1 and  from  level  3 to  level  1. 

In  other  specific  cases  which  involve  stepwise  fluorescence  processes, 
it  is  possible  to  make  different  assumptions  which  can  be  used  to  simplify 
the  form  of  hq.  (45).  However,  because  stepwise  fluorescence  is  not  yet 
widely  used  in  analytical  atomic  spectrometry,  these  cases  were  not  indiv- 
idually treated. 
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l UNt:i  UNIONS 

The  model  for  saturation  of  atomic  levels  presented  in  the  preceding 
section  is  far  more  general  than  those  ordinarily  applied  to  saturated  AFS 
(.1-10,  18,  19).  Not  only  does  the  model  lend  itself  to  the  convenient  des- 
cription of  saturation  effects  in  non-resonance  fluorescence,  it  also  is 
applicable  to  atomic  systems  driven  toward  saturation  by  transient  pulses. 

With  this  theoretical  treatment,  it  becomes  possible  to  predict  the 
effect  of  various  experimental  variables  on  the  saturation  process.  As  an 
example.  Figure  8 displays  the  predicted  dependence  of  saturation  spectral 
power  density  (pg)  on  the  duration  (tj)  of  the  exciting  radiation  pulse. 

The  case  shown  is  the  Ba  transition  at  553.5  nm,  with  an  assumed  lifetime 
(t)  of  8.5  ns  and  a collisions!  quenching  rate  equal  to  lf)A2i.  Consequently, 
the  true  excited  state  lifetime  ( t ' = 1/(A.  i + -?]))  is  approximately  0.8 
ns.  One  can  see  from  Figure  8 t hat  the  saturation  spectral  power  density 
increases  markedly  as  the  excitation  pulse  becomes  narrower  than  the  true 
lifetime,  but  becomes  nearly  constant  for  excitation  pulses  much  longer  than 
the  lifetime. 

No  attempt  has  been  made  to  exhaustively  illustrate  the  applicability 
of  the  model  to  non-resonance  AFS.  Rather,  the  two  most  commonly  employed 
non -resonance  methods  (direct-line  and  stepwise  fluorescence)  were  selected 
to  show  how  pertinent  excited-state  populations  could  be  predicted  and  how 
these  populations  relate  to  exciting  source  powers  necessary  to  reach  satur- 
at ion . 

The  calculation  of  source  power  densities  necessary  to  produce  satura- 
tion is  perhaps  the  most  significant  capability  of  the  proposed  model.  With 
such  a capability,  it  should  be  possible  to  predict  the  spectral  power  den- 
sity which  would  generate  saturation  effects  of  the  kind  found  to  be  so  ben- 
eficial to  analytical  AFS.  As  we  hope  to  show  in  a later  communication, 
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accurate  measurement  of  these  required  spectral  power  densities  is  extremely 
tedious  and  must  be  performed  with  great  care.  With  knowledge  of  these 
requisite  power  levels,  it  should  be  possible  to  assess  the  feasibility  of 
approaching  atomic  saturation  with  conventional  (non-laser  sources. 

As  with  all  theoretical  models,  the  one  proposed  herein  must  be  supported 
by  experimental  evidence.  A study  to  provide  this  evidence  has  already  been 
completed  and  will  be  reported  shortly.  In  that  study,  the  proposed  model 
has  been  found  to  agree  far  better  than  earlier  models  with  experimental  re- 
sults and  can  be  used  to  determine  values  and  trends  relating  to  the  satura- 
tion of  atomic  energy  levels. 
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a.  = \y  i *•  i = 1/ i 1 for  2- level  atomic  model 

b?  = H i ^ i for  2-level  atomic  model 

3e  »b  )e»o  <e  - lumped  parameters  referring  to  3-energy  level  model  with  the 
assumption  of  thermal  equilibrium  between  levels  3 and  1. 

See  l.qs.  36-38. 

1 isil’  si1  is  = lumped  parameters  pertaining  to  direct-line  fluorescence  from 
a 3- level  atomic  model  in  which  the  lower  level  involved  in 
the  transition  is  assumed  to  have  a steady-state  population. 

See.  l.qs.  4 2-44. 

anr»bnr  ~ lumped  parameters  pertaining  to  stepwise  fluorescence  from  a 3- 
level  atomic  model  in  which  a steady-state  population  is  assumed 
for  the  upper  level  in  the  fluorescence  transition.  See.  F.qs.  48-49. 

Aj i = Einstein  coefficient  for  spontaneous  emission  from  level  j to  level  i 
( sec"  1 ) 

Bp  = detected  fluorescence  radiance 
c = speed  of  light  (3  x lu!0cm/sec) 

D = departure  of  rate-controlled  average  population  of  excited  level  (under 
non- steady-state  irradiation]  from  the  population  produced  by  steady- 
state  irradiation.  See  Lq.  18. 

hi  = Energy  of  atomic  level  i 

f-  = Fractional  (as  %)  error  in  saturation  spectral  power  density  caused  by 
erroneous  assumption  of  steady-state  irradiation. 

g^  = degeneracy  (statistical  weight)  of  the  ith  energy  level 

h = Planck's  constant  = 6.624  x io-  i4  J-sec 

k = Boltzmann  constant 

A j i = wavelength  corresponding  to  a transition  from  energy  level  j to  energy 
level  i 

ni(t)  = Number  of  atoms  in  energy  level  i at  time  t 

N = total  number  of  atoms  in  aJl J_  levels  under  consideration 

= Steady-state  population  of  upper  level  2 produced  by  continuous  ir- 
radiation 

Nj  = nj (t)  = average  number  of  atoms  in  upper  level  j during  the  excitation 
pu  1 se 

(Nj)sat  = population  of  level  j under  conditions  of  complete  saturation. 

vji  = spectral  frequency  corresponding  to  a transition  from  level  j to 'level  i 

p 3es  = saturation  spectral  power  density  corresponding  to  direct-line  fluro- 
rescence  from  a 3-energy  level  model  atom  with  assumed  thermal  equil- 
ibrium between  levels  1 and  3 


p(t)  = spectral  power  density  of  incident  beam  (watts/cm^Hz) 
(io  = peak  spectral  power  density  of  exciting  radiation 

Ps 


a 


j-i.  = saturation  spectral  power  density  for  2-level  atomic  model 

,sNSS  = peak  spectral  power  density  of  a transient  (non-steady- state)  exci- 
tation pulse  which  produces  an  excited  level  population  half  that 
generated  during  saturation.  Termed  the  non-steady-state  saturation 
spectral  power  density. 


t = time 

to  = excitation  pulse  length. 

T = absolute  temperature  (°K) 

t ' = observed  fluorescence  lifetime  for  transition  of  interest  = l/a? 
Z = M3/a3i  (see  Eq.  34) 


Zij  = rate  of  collisional  (radiationless)  transition  from  level  i to  level 
i (sec-1) 


Z2l/A21  = quenching  factor  for  a 2-level  model 
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Table  m.  Dependence  of  calculated  saturation  spectral 
power  density  on  the  assumed  value  of  the 
fluorescence  lifetime 

Assumed  Lifetime  of 
Atomic  Transition 
T (ns) 

k 

e 

1 6 

a.  The  parameter  e expresses  the  per  cent  error  in  the  predicted 
saturation  power  density  which  would  be  caused  by  assuming  a value  for 
the  fluorescence  lifetime  (t  ) that  is  50$  in  error. 


Table  IV  Values  of  saturation  power  density  predicted  by  the 

two  energy- level  model 
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Figure  1.  Schematic  diagram  for  excitation  and  saturation  of  a two  energy 
level  atomic  system. 

A)  Energy  level  diagram  bowing  transition  rates;  Bj2  and  B2  j 

are  the  Einstein  coefficients  for  induced  transitions  (1  ■+•  2 
and  2 - 1)  respectively,  A2 i is  the  Einstein  coefficient  for 
spontaneous  emission  (2  ■*  1),  Z2 i is  the  rate  for  the  colli- 
sionally  induced  transition  (2  1)  and  p(t)  is  the  spectral 

power  density  of  the  exciting  beam  of  radiation,  a2  = A2 l + Z2 j 

B)  Boxcar  function  assumed  to  describe  the  time  behavior  of  the 
exciting  radiation;  to  and  po  are  the  width  and  the  peak 
spectral  power  density  of  the  radiation  pulse,  respectively. 

Figure  2.  Time  behavior  of  the  fraction  of  excited  atoms  (solid  curve) 

during  application  of  a 5 ns  duration  pulse  of  exciting  radiation. 
Pertinent  assumed  rate  constants  governing  excitation  and  decay 
are  given  in  the  text. 

N?  _ fractional  population  of  excited  state.  Rise  and  fall  times 
N calculated  at  1/e  points. 

Figure  3.  Fraction  of  excited  atoms  produced  by  non-steady  state  irradia- 
tion by  a laser  of  varying  spectral  power  density.  Degrees  of 
excitation  for  various  pulse  lengths  (to)  are  shown. 

A)  to  = 1 us 

B)  t o = 10  ns 

C)  tg  = 1 ns 

Fraction  of  atoms  in  excited  state  when  laser  power  density  is 
infinite  = 0.75  (see  text  for  explanation). 

Figure  4.  Predicted  saturation  curves  showing  the  effect  of  quenching  col- 
lisions on  excited-state  population.  Excitation  pulse  is  of  1 
ps  duration.  Curves  A,  B and  C correspond  to  atom  environments 
in  which  the  rates  of  collisions  are  0,  5 and  10  times  the  rate 
of  spontaneous  emission,  respectively.  N2/N  is  the  fraction  of 
all  atoms  in  the  upper  state;  for  complete  saturation,  N2/N  = 
0.75. 

Figure  5.  Predicted  saturation  curves  showing  the  effect  of  quenching  col- 
lisions on  excited-state  population.  Excitation  pulse  of  6.5  ns 
duration.  Curves  A,  B,  and  C correspond  to  atom-environments  in 
which  the  rates  of  collisions  are  0,  5 and  10  times  the  rate  of 
spontaneous  emission,  respectively.  Compare  to  Figure  4. 

Figure  6.  Effect  of  the  accuracy  of  the  values  of  fluorescence  lifetime  on 
the  saturation  curves. 


A)  x = 4 ns 

B)  t = 8 ns 

C)  t ■■  16  ns 


Figure  7.  Schematic  model  for  excitation  and  saturation  of  a three-energy 
level  atomic  system. 

A)  Energy  level  diagram  showing  the  activation  and  deactivation 
pathways. 

B)  Boxcar  function  assumed  to  describe  the  time  behavior  of  the 
exciting  radiation. 

See  text  for  discussion  of  rates  of  various  transitions. 


Figure  8.  The  predicted  effect  of  excitation  pulse  length  on  saturation 
spectral  power  density  for  the  Ba  553.5  nm  transition.  True 
lifetime  of  the  transition  assumed  to  be  0.8  ns. 
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